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Abstract. - We show that Kramers rate theory gives a straightforward, accurate estimate 
of the closing time r c of a semiflexible polymer that is valid in cases of physical interest. The 
calculation also reveals how the time scales of chain relaxation and closing are intertwined, 
illuminating an apparent conflict between two ways of calculating r c in the flexible limit. 



The looping of polymers is a physical process that allows contact and chemical reaction 
between chain segments that would otherwise be too distant to interact. Polymer loops are 
particularly important in biology: In gene regulation, looping allows a DNA-bound protein to 
interact with a distant target site on the DNA, greatly multiplying enzyme reaction rates 
Similarly, DNA looping in the 30-nm chromatin fiber may trigger the initiation of DNA 
replication at different sites along the DNA by enabling long-distance interactions 3 . In 
protein folding, two distant residues start to come into contact via looping 1 1I5| . Measurements 
of loop formation in single-stranded DNA segments with complementary ends have also been 
used to extract elasticity information (e.g., the sequence-dependent stiffness of single-stranded 
DNA 0). 

Despite its importance and despite considerable theoretical effort, there are relatively few 
analytical results concerning the dynamics of loop formation. Even for the simplest case of 
an ideally flexible polymer with no hydrodynamic effects (or simply a Rouse chain), there 
are two rival theoretical approaches that lead to contradictory results: Szabo, Schulten, and 
Schulten (SSS) conclude that the time for a loop to form ("closing time" r c ) should scale for 
moderately large polymer lengths L as r c ~ L 3 ^ 2 J?|, while Doi, applying Wilemski-Fixmann 
(WF) theory finds Tuoi ~ L 2 [3]. The discrepancy between the two continues to spur 
debate |1UII11| . For the important case of stiff chains |12II13| . where the polymer length L is 
not too much longer than the persistence length £ p , only limited numerical results are known 
[see, for example, |14I15| and references therein]. The main difficulty arises from the interplay 
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between two seemingly distinct processes: chain relaxation and chain closure. This interplay 
is unique to a polymeric system and originates from the chain connectivity of a polymer 
immersed in a noisy environment. 

In this article, we argue that Kramers' rate theory i6.,,17j applies to the most physically 
relevant cases and leads to analytical results for r c . We capture, for the first time, that 
there is a minimum loop- formation time for chain lengths of approximately 3-4 £ p . Roughly 
speaking, shorter chains require too much energy relative to the thermal energy fc^T, while 
longer chains need to search too many conformations for ends to "find" each other. We also 
show that consideration of the requirements for Kramers theory to apply leads one naturally 
to identify different regimes governing the closing time r c . This classification shows how the 
physics of chain relaxation is intertwined with that of chain closing and clarifies the above- 
mentioned controversy between the SSS and Doi approaches to loop-formation dynamics. 

Consider a chain of length I = L/£ p with two ends that react when first brought within 
a distance a of each other ( "diffusion-limited" loop formation dynamics) . We apply Kramers 
rate theory, viewing the process as a noise- assisted "tunneling" over a potential barrier. Af- 
ter first presenting the straightforward calculation, we then consider carefully its domain of 
applicability, followed by a scaling description of loop formation outside this domain. 

The basic idea is to project the internal degrees of freedom of the polymer chain onto 
a single "reaction coordinate" r = R/£ p , with R the end-to-end distance of the chain. The 
reduced, one-dimensional dynamics then obey a Langevin equation of the form 

§ = -^W',0-K(t), (1) 

where D = 2Dq is twice the diffusion coefficient of a monomer (both ends diffuse JU]) and £(t) 
represents Gaussian white noise: (£(t)) = and = 2D5{t — t'), with (...) a thermal 

average. The dynamics are governed by an effective potential U(r,£). Strictly speaking, this 
description is valid for a chain in local equilibrium, for which we can write 

U(r,£) = -k B TlnP{r,£), (2) 

where P(r,£) = 4irr 2 G(r,£) is the radial distribution function of end-to-end distances r of a 
polymer of length £ and G(r,£) = GQre — r \;£), the angle-averaged distribution function for 
the end-to-end vector r = Tt — Vq. We assume isotropic chemical interactions between end 
monomers, so that end binding can be modeled by adding to U a smooth short-range potential 
f(r/a), with a = a/£ p the interaction range. A typical distribution function and resulting 
effective potentials are shown in fig. Q 

Because polymers - whatever their stiffness - have a most probable end-to-end separation 
(radius of gyration) , there is a local minimum in the effective potential at r?, (bottom) , which is 
~ £ in the stiff-chain limit and ~ ^fl in the flexible-chain limit, neglecting self-avoidance effects. 
Also notice the barrier to chain closing at r t ~ a (top), which is created by the balance of 
chain entropy and bending energy as implied by U (r, £) . The short-range attractive potential 
then rounds off the barrier. The resulting effective potential has thus the qualitative form 
often assumed in Kramers-rate calculations. 

In the limit of strong damping |18j . the time needed to tunnel over the barrier (mean 
first-passage time), calculated using Kramers rate theory, is 



TKr = 



D ujtLUh ( AU 



(3) 
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Fig. 1 - (a) The radial distribution density P(r,£ = 3). The dashed line shows the effect of a short- 
range interaction between the two polymer ends, (b) The resulting effective potential of the chain. 
Arrows denote the top and bottom of the effective potential well, as used in the Kramers calculation. 



where the well curvatures ui(r) = j-y/ d rr ll(r, £) are evaluated at the top and bottom of the 
effective potential U(r,l). The exponential factor is 

/ AU\ P(r u l) a 2 G(Q,£) , 
We find the surprisingly simple result |19) 

TKr{i)=C^^ (5) 



aD G (t)' 



1 I 2 

a dimension- 



with Go = G(0,£) and with C (r b ,G(r b ,£)) = V2irr 2 b G(r b , £)/ - 

less prefactor that is practically a constant for all £ (see below). 

Eq. [5] shows that the closing time may be estimated using the static distribution G(r,£). 
Unfortunately, no analytic expression for G(r, £) has been found that is accurate for all r 
and £, and one must make a pastiche of approximations. For r = and i < 10, we use an 
approximation for a wormlike chain derived by Shimada and Yamakawa |2()U21j : Gq(£) = 
(896.32/£ 5 )exp(-14.054/^ + 0.246^). Note here that the l/£-term in the exponent solely 
arises from bending energy, while the rest comes from chain fluctuations around the lowest- 
energy conformation. For r = and large £, we use an interpolative formula due to Ringrose 
et al. that blends SY with the result for a freely jointed chain, Gq{£) ~ £~ 3 / 2 Near 
r = r b , we use an approximation derived by Thirumalai and Ha (TH) |23| . valid to 10%: 
G(r,£) = n(£)[l — (|) 2 ]~2 exp { — T [1-^r/ty 1 } }' tnc normalization factor n{£) fixed by 

requiring j 47rr 2 G(r, £)dr — 1. Note that a more accurate but more complicated expression 
recently derived by Winkler j^] gives essentially the same results. Using TH, we find that the 
dimensionless prefactor C(£) of eq.[5]is O(10 _1 ), varying less than a factor of 2 over < £ < 00. 

In fig. we plot the TKr W that results from eq. [S] using the various approximations 
to G(r, £) discussed above. The solid curve uses the Ringrose expression for all i, while the 
dashed curved uses SY for small £. The two curves compare well with recent simulations using 
parameters appropriate to double-stranded DNA ^3]. Note that the material parameters of 
the simulation were used (see caption). Considering the heuristic nature of the arguments, 
the agreement is excellent. 

One striking feature of the plot of T^ r (£) is the existence of a minimum at £ as 3.4, where 

^ = 0-78 A. (6) 
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Fig. 2 - Closing time r c vs. chain length, (a) BD simulation |14| (empty circles) and Kramers theory 
(eq. [SJ are shown. For direct comparison, we used the same parameters as in the Ref. |14| (bead 
size = 3.18 nm for D = 2D = 1.54 x 10" n m 2 /s and a = 0.1) with £ p = 50nm. For Go(£), we 
used the SY result 12(11 and an interpolation [221 (see text). Relaxation times tr for these parameters 
are also shown (triangular symbols), with the £ 4 and £ 2 scaling regimes apparent in the inset, (b) 
Single- "particle" MC simulations of t c with the potential U/ksT = — log[P(r, £)] taken from Fig.^t. 
Here, r c is a first-contact time averaged over about 2000 realizations of the initial position randomly 
selected from P(r,£). We have chosen a = 0.25, 0.5, 0.75, 1.0. As expected, r c ~ (inset). 

In eq. the prefactor 0.78 is calculated by Monte Carlo simulation of G(r,£), in units of 
seconds, and is about 10% less than the prefactor obtained using the TH approximation. As 
mentioned above, the existence of a minimum in tkt reflects a balance between the energy of 
bending and the entropy of conformations that must be searched for two ends to meet. 

For the above Kramers-rate calculation to hold so that Tjc r equals r c , three conditions 
must be satisfied: The damping must be sufficiently strong; the barrier height AU must be 
large compared to ksT; and the global chain relaxation time tr, a characteristic time scale 
for chain deformation, must be much shorter than the Kramers time tkt- 

The first condition is normally satisfied for molecules in solution |18j . For the second, since 
there is a minimum in the effective potential at r^, we require that a be sufficiently less than 
Tb so that the barrier height is large. The condition AU /fcsT = 1 is shown in fig.[3]as a dotted 
line in the £ — a parameter plane, using a diffusion constant appropriate to double-stranded 
DNA. To the left of the dashed line, the barrier height is larger than fc^T. 

The third condition, tr <C t#>, is more subtle and requires discussion. In using a "one- 
particle" description of chain closing dynamics, we are assuming that all internal degrees of 
freedom of the polymer chain have relaxed. As a result, the end-to-end distance is the only 
dynamic variable (Cf. eq. ^) . This assumption of local equilibrium allows one to apply the 
equilibrium distribution function G(r, £) and implies that the effective potential derived from 
G is time independent. If the chain relaxation times are too long, the potential effectively 
becomes time dependent and has to be obtained self-consistently, along with the motion of 
the internal modes. We thus compare the scaling behavior of tr{£) with TKr{() and t c (£) in 
both the flexible (£^> 1) and stiff-chain {£ < 1) limits. 

In the flexible limit, we can use the Rouse model to estimate the longest relaxation time, 
which gives tr ~ £ 2 , in units of the basic time scale ip/D. By contrast, at large £ eq. Ogives 
i~Kr ~ £ 3 ^ 2 /a. (This is just the result of SSS |7I10| and has been confirmed by single- "particle" 
simulations — see fig. |2Ib) and the caption.) Thus, when £ > 1/a 2 , the third condition is 
violated and the Kramers calculation does not hold. Nonetheless, we can still estimate the 
upper-limit of r c : The closing time is at most a time necessary for the slowest "random 
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Fig. 3 - Scaling regimes in the £-a for DNA (see text). Region I is the Kramers regime, with r c > tr; 
Region II is the dynamic-fluctuation regime. In the primed regions to the right of the dashed line, 
AU/ksT < 1. The black region is unphysical: a> L. 



"walker" to travel, with diffusion constant that of the entire chain D c hain ~ D/l, the end-to- 
end distance r. Since r ~ we have r c < r 2 /D c hain ~ £ 2 /D ~ tr. In other words, when 
the third condition does not hold, t c is not tk> but is set by the Rouse time tr. 

In the stiff limit, the physics is dominated by the bending energy of a rod |25) . leading 
simultaneously to faster relaxation times tr and higher energy barriers, which implies that 
the Kramers calculation should be valid. To see this, recall that the bending energy of an 
elastic rod is, by symmetry, proportional to the square of the rod curvature. Since the lowest 
energy corresponds to a uniform curvature of radius 1Z, the bending energy near the rod limit 
E b /k B T = \t p L/Tl 2 ~ 2£ p (L ~ R)/L 2 . 

Thus, if we track the relative separation R of the endpoints of a thermally excited rod, it 
behaves like a particle subject to a constant restoring force f c = 2fcsT£ p /i 2 . The appropriate 
Langevin equation for R is then of the form R+ D ^ T n f c S = £ c /j j n (i), with £ c hain(t) the 
random force. This implies that the time to relax a distance of order L is ksTL/ D c hainfc = 
I? /2£ p D chain- Since the rod moves coherently, the diffusion coefficient of the chain Dchain ~ 
Db/L, leading to tr ~ ^TbB ' wnere & i s the monomer size. As a result, for L < l p , the third 
condition (tr <§; Tkt) is always satisfied: the lower-limit of r c is given by a time scale for a 
random walk to travel a distance R^L, thus TKr~T c > R 2 / Dchain^L? /bD > L 4 / £ p bD>TR. 

To summarize, tr ~ L A for £ < 1 and ~ L 2 for £ ^> 1. Thus, for large enough i, tr 
becomes larger than the Kramers estimate as shown in fig. |2K and in the inset. In fig- 01 
we plot t r (£) — r Kr (l) in the £— a plane. The white area is Region I (Kramers Regime), where 
i~Kr > tr, and therefore r c ~ tkt- The shaded area is Region II ("dynamical fluctuation" 
Regime, see below), where tkt < tr, and therefore t c ~ tr. Areas I' and II' show where 
At/ < UrT. The black region, defined by a > £, is unphysical. 

In Region II, the relaxation and closing processes are coupled. In this case, one may have 
to solve an A^-particle diffusion problem, subject to a boundary condition that is difficult to 
impose [HJEIEII- Nevertheless, much insight can still be obtained from the simple scaling 
analysis of random walks given above. In this view, a chain can close because the two ends 
randomly meet each other while freely relaxing. The existence of such a regime, where r c ~ tr, 
is a unique feature of flexible chains (fig. that we denote the "dynamical fluctuation" 
regime — the dynamic fluctuation 6R(t) = v/ ([R{t) — R(0)] 2 ) grows up to R as t — > tr and 
thus can assist chain closing. For a Rouse chain, 5R(t) can be given as a sum of Rouse 
modes |28] and, in our simple scaling analysis, r c can be inferred by analyzing this. Short- 
time behavior of 5R(t) reflects the internal motion and varies as SR(t) - y/i for t < tr. We, 
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however, argue that this will not appreciably influence r c , as SR(t) — > R only when t — > tr. 
In other words, r c is governed by the slowest mode and our assertion of r c ~ tr will not 
be invalidated by the internal motion, which is important at time scales much smaller than 
t c (or tr). In the stiff-chain limit, this dynamical fluctuation regime disappears. Note that 
the boundaries between Regions I and II are not sharp but are crossovers. Loop-formation 
kinetics in the crossover area will likely combine aspects of both regimes, as indicated in recent 
simulations jl(J| and by results that show that tsss and tdoi are respectively lower and upper 
bounds for r c Similarly, based on their BD simulation results, Podtelezhnikov et al. |27| 
suggested that r c ~ tr/cx near the boundaries. 

Our discussion has neglected hydrodynamic effects and excluded- volume interactions. Both 
can influence chain relaxation and closing simultaneously. The hydrodynamic effect will not 
change TKri since it is a function of the equilibrium distribution G(r, I). However, the hydro- 
dynamic interaction tends to promote chain relaxation (e.g. in the Zimm model, tr ~ £ 3 / 2 , in 
contrast to tr ~ £ 2 in the Rouse model considered here |2E|) by increasing the mobility of the 
chain, resulting in a wider Kramers regime than implied by Fig. |31 On the other hand, the 
excluded- volume interaction both decreases D 'chain and reduces Go |28H29j . But for loops of 
just a few persistence lengths, which are the most physically relevant (see below), both effects 
are expected to be minor. A final caveat is that we have assumed isotropic binding inter- 
actions. While mathematically simpler and relevant to simulations |14|. most real polymers 
have directional bonding. In the Kramers calculation, this would modify Gq(£). 

The Kramers calculation holds in Region I of the £ — a parameter space shown in fig. [3] 
What are the physically relevant values of a and 11 The interaction distance a — a£ p will be 
the thickness of the polymer, or less. For polymers of biological interest, the persistence length 
will be typically at least this size and often much larger. For example, for double-stranded 
DNA, the monomer size is 0.34nm while the persistence length is 50 nm. For chromatin, the 
thickness is 30 nm, comparable to its persistence length [^0]. Thus, we generally expect a < I 
and sometimes a«l. 

What are the relevant values of 11 Although polymers in principle may have any length, the 
existence of a minimum closing time T^ r (eq. |HJl leads one to speculate that where looping is 
biologically relevant, polymer lengths near £ ss 3—4 might be favored because they minimize r c . 
In this regime, the Kramers calculation will be valid, for small a. Thus, biological selectivity 
may arise from a physical mechanism. For example, a recent study of Jun et al. 3 on DNA 
replication noted that the typical spacings between replication origins in early embryo Xenopus 
are 3-4 times the £ p of chromatin, the DNA-protein complex present during replication. It 
is then natural to speculate that origins are related by looping and that the spacing may be 
selected to maximize the contact rate of origins, optimizing replication efficiency. 

In conclusion, we have shown that Kramers rate theory gives a straightforward estimate 
of the closing time of a semiflexible polymer. Although phenomenological, the calculation 
explains the existence of a minimum closing time and accurately reproduces numerical sim- 
ulations. Moreover, considering the requirements for the calculation to hold shows how the 
intertwining of the relaxation time with the closing time explains the apparently conflicting 
results for r c (SSS and Doi). Fortunately, the physically relevant cases are precisely the ones 
where the Kramers calculation is expected to hold and may even be selected biologically 
through evolution. 
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